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Abstract 

We present an exact solution of headway distribution of the asym- 
metric simple exclusion model with open boundary conditions and 
compare it to the headway distributions of the highway traffic. 

1 Introduction 

Recently, a new tendency has appeared in the physics of traffic systems. Besides macroscopic 
characteristics of traffic (flow, density, etc.), researchers are focused on microscopic properties 
as well (0). One of these properties is the distribution of the distances (or time intervals) 
between neighboring cars, so called headway distribution. This function was investigated for 
existing traffic models and compared to the various regimes of the real traffic. It was observed 
that the real traffic data as well as many traffic models lead to a headway distribution which 
resembles the spacing distribution of the random-matrices eigenvalues (@,|§],|H|). To verify this 
observation, we investigate the asymmetric simple exclusion process (ASEP) which is the well 
known exactly solvable model of traffic (@,|§,||6|). We find the exact formula for headway 
distribution of the asymmetric simple exclusion model and compare it with the eigenvalues 
distributions obtained from random matrix theory. 

In the traffic theory, the spacing distribution which was often called headway distribution 
P(s) is defined as probability distribution of the distances s between neighboring cars (particles 
of the model). For the sake of simplicity we assume that (s) = / °° s P(s) ds = 1. Among many 
others, we discern two special traffic systems in real roads. The first is the bus transport where 
no fixed schedule exists (see @],|§J) and where the interaction between bus drivers exists owing 
to the check points in which the information about the time interval to the previous bus is 
transmitted to a driver. It was found that the headway distribution of that bus transport 
behaves as 

P(s) =As 2 e- Bs \ (1) 

where A,B are constants determined by conditions f£°P(s)ds = 1 and s P(s) ds = 1. 
The equation (||) originates from random matrix theory where it describes the distribution of 
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distances between eigenvalues of the unitary random matrices. 

The second system described here is highway traffic for which we distinguish three basic regimes 
(§): 

• a) free-flow regime, where the density of cars (number of cars per unit 
of length) is low and the cars are therefore moving as free. Headway 
of this regime is reduced to the Poissonian distribution 

P (s)=e- S (2) 

• b) synchronized regime, where the mean velocity is considerably re- 
duced compared to the free-flow states due to a higher density of cars, 
but cars are still moving (no traffic jam). Headway distribution of 
these states is well described by the formula 

p ^ ) = {j W^Y sUe ~ {v+1)s ^ (3) 

where v is a free parameter and T is Gamma function. 

• c) stop-and-go regime, where cars are moving in traffic jams 

Note that it is possible to obtain the formula (|2|) taking the limit v — > + in formula (|||). 
Distribution P\(s) used to be called as semi-Poissonian statistics. 

Roughly, particles in the free-flow regime do not interact between each other, contrary to 
the synchronized states when the interaction leads to the particle repulsion and to the fact that 
lim s ^ + P u (s) = for positive v. One of aims of this article is to answer whether it is possible 
to get a synchronized state for the asymmetric simple exclusion model by an appropriate choice 
of input parameters, or the model is in the free regime for all parameters. 



2 ASEP with open boundaries 

Consider sites located on a chain of length N, Each site i (1 < i < N) is either occupied by 
a particle, or empty. During the time interval dt, each particle jumps to the next site (in the 
defined direction) with probability dt if the target site is empty and does not hop otherwise. 
Moreover, during the time interval dt, a particle can enter the first chain position (i = 1) with 
probability adt (a £ (0, 1)) if this position is empty and a particle at the site i = N can leave 
the chain with probability f3dt (J3 € (0, 1)) if this site is occupied. These rules define a simple 
traffic model with a hard-core repulsion which decreases the probability for occurrence of the 
two particles close each other. This was de facto the reason why we have chosen ASEP for our 
investigation. 

Two alternative methods of exact solution of the model are outlined in [|]]. We will use the 
matrix method in our computations. Consider the matrices D,E and vectors (w\ and \v) 
satisfying the following algebraic rules 

DE = D + E (4) 
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(w\E = -(w\ (5) 
a 

D\v) = hv) (6) 

Then we can determine the probability of the arbitrary configuration C = (ni, mi, n,2, m.2, • • •) 
of the model through the formula 

1 (w\D ni E mi D n2 E m2 ... \v) 
Z N (w\v) 

where n\ is the number of the empty sites, starting with the entering site of the chain (and 
(ni + l)-th site is occupied), mi is the number of occupied sites starting with the site number 
ni + 1 (and (ri\ +mi + l)-th site is empty), and so on. Zn is a normalization constant. Hence 
J2i n i + J2j m j = N. Owing to the fact 



m (In — m — 1)1 ,■ , 

(D + E) n =Y\ —^T, rr - ^ V E l D r 

v ; ^ n\(n-m)\ f-^ 

m=l v y i=0 



we can rewrite the normalization factor Zn to the form 



in— i 



_ (w\(D + E) N \v) _ " m(2N-m-l)\ ^fiy ( 1 \ 

(w\v) ^ N\(N-m)\ ^\a) \p) W 



Hence, 



_ m (2N - m - 1)1 p-™' 1 - a'™- 1 
Zn ~ \ N\(N-m)\ - a-i (10) 

m=l y i ' 

These are the known results 01 which will be used below. 



3 Headway distribution in the ASEP (a special case) 

Our aim is an analytical expression for the spacing distribution, which means the probability 
of the gap of the length k behind the particle being on the i-th site. One can determine this 
probability PW (k) from the equation (0) as follows 

pW(jfe) = —j— { w \{D + E) i - 1 DE k - 1 D{D + E) N - i - k \v), (11) 
Zn \w\v) 

since the first i — 1 sites could be arbitrary occupied or not, i-th site is occupied, the following 
k — 1 sites are empty, the site i + k is occupied and the rest sites could be arbitrary occupied 
or not. If these rules are satisfied, one can see the gap of the length k behind the i-th position. 
The constant A can be calculated from the condition 

N-i 

E p(i) (*o = i 
k=i 

Generally, the matrices D,E are infinite-dimensional, but provided the condition 

a + (3=1 (12) 
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one can find D, E to be real numbers. The equations (||) and (]6|) imply that D = 1/(3, E = 1/a 
and the rule (||) imposes that 1 / a + 1 / / g = 1 / a p, which leads to the equation dl2j). Instituting 
(0) and @ into the @ we obtain 

=^a 2 /3 fc ~ 1 (A: > 1) (13) 

Strictly speaking, the derivation of the formula (^) from the equation is possible 
only for fc > 3 since for k = 1,2 the relation (|Tl|) has a different form. But it can easily be 
verified that (|i~3|) holds true also for A; = 1,2. Note now that the result (|i~3|) is independent 
of the value i For large N, after the re-scaling to (s) = 1, equation (|T^ ) changes to the form 
-P(s) := P(k — 1) = e -5 , which is the Poissonian distribution. This clearly shows that under the 
condition (|i~2|), the interaction among the particles vanishes in the continuous limit and hence 
the movement of the particles is asymptotically free (in the sense mentioned in the section 2). 
Also the real headway distribution (without the re-scaling) can be determined in this section. 



If (pi) holds true then 

N-2 / i \ N+l-k 



/ i \ / 1 \ iV + l-« 



,13 

when the normalization pre-factor M can be computed either from the condition J2k=i (^) = 
1 (through all the lengths of the gap) or from the equation 

1 _ (w\(D + E) N - i D{D + E) N ~ i \v) - (w\(D + E) i - 1 DE N - l \v) 

which takes into account all the possible combinations with occupied site number i and rejects 
the combinations when all the sites (behind the site number i) are empty. This provides 

1 / 1 \ N ~ i ( 1 x N 



which together with (|14| ) yields 



In the large N limit, one can calculate from ( |T2"| ) that 

P {i) {k) = a(3 k ~ 1 

and hence 

oo oo 
k=l k=l 

for the mean headway. 



4 Headway distribution of the ASEP (general case) 



In this section we are solving the ASEP (equation (]TT|) , strictly speaking) for arbitrary values 
of the parameters a, (3. Introduce now the following pre-factors 
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uj = ui(m, k, i, N,p) 



mp(2N -2i-2k-m- l)!(2i - p - 3)! 
(N — i — k)\(N -i-k- m)\{i - - p - 1)! 

A 1 = A 1 (p,q,a,w) = I y p - q -i ) 



A 2 = A 2 (p,q,b,w) 



p — q — b + w — 1 
w — 1 



With the help of them and the computations presented in an Appendix we obtain the 
desired result 
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Re-scaled headway s 



where 



FIGURE 1. Headway distribution of the ASEP . The curve 
is (^). Bars and plus signs show the result ( |l5|) after the re- 
scaling to mean headway equal to one for a = 0.15, /3 = 0.6 
and a = 0.3, fi = 0.5, respectively. 



A 



i— 1 N—i—k m 



P {i) (k) = ^ L Y. E E u(S 1 + S 2 + S 3 + S± + S 5 ) (15) 

An p =l m =l 1=0 



Si 



p-1 



p— 1 fe— 1 / to 
g=0u>=l \tt=l 



^ \ q+a / i \ m— l+l PI / 1 \ (J /^\m+f)-i+l N 



/3 



p— 1 2 + fc — 1 / UJ 



'E E 

q=0 ui=l \ a=l 



i y +a / i 

a) \J3 



m-l V~1 
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*-(=) G) + £G) ® 

/ 1 \ P / 1 \ W— i+l z+k—l / -1 \ p+uj / i \ m— Z 

G) +i £G) G) 

The equation (|15[) presents the final result of this article. In the special case ( O ) the equa- 
tions (HH) and (|lj) provide the same values, which supports the correctness of our result. Now, 
one can show that the headway distribution of the ASEP (after re-scaling) is of a Poissonian 
type for any parameters a, (3 (see Figure 1). 

5 Conclusion 

To conclude, the ASEP shows the Poissonian behavior for all possible values of the parameters 
a and (3. Thus, this model can simulate only the traffic in free-flow regime, which means that 

lim P(s) » (16) 

To be precise, we have to express that as the parameter [3 goes to zero ((3 — » + ) the Poisso- 
nian distribution changes since one can easily to follow that P(s) = S(s) = 5(k — 1) for (3 = 
and This holds true due to the shift of the probability P{s) := P(k — 1). The above 

mentioned shift is done because P(k) = for every k < 1 in the original model. However, for 
every [3 we get the headway distribution satisfying the condition (|l6|). 

On contrary, in a different variant of the ASE-model (high speed traffic model with open 
boundaries) is presented. The particles of this model also enter the system at the left and leave 
at the right of the chain. Moreover, all particles will move with their maximal possible speed, 
which is given by the speed limit v max (a maximal length of the jump). For v max = 1 one 
can obtain the model investigated in this article. Numerical simulations of the high speed 
model with v max > 2 clearly show that the gap distribution can start from the origin (i.e. 
lim s _ >0 + P(s) = 0) for some values a, (3. 

Hence, the high speed model has the synchronized states as well, on contrary to the asymmetric 
simple exclusion process and it can be used to the full simulation of the real-road traffic. 
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7 Appendix 

At the beginning, we cite some useful lemmas. 
Lemma 1 

Let m, n, a are the natural numbers and a < m. Then 

Y^fn — i + rn — l\_fn + m — a\ , , 

Proof: by means of the mathematical induction for the numbers m and n 



Lemma 2 

Assume that D, E are arbitrary matrices fulfilling the equation DE = D + E. Let m, n are 
the natural numbers. Then 

D"ET=t{ n ~ " 1 ) * + t ( " " " 1 ) & US) 

Proof: by means of the mathematical induction for the numbers m and n and ( |l~7| ) 
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Now we can rewrite the equation (11) in order to use the rules (H),©. Thus, it is necessary 
to replace the matrices D on the right-hand side and matrices E to the left-hand side of the 
equation (fll]). The expression (@) yields 



(D + E) 



v _._ fe N z^ k m (2N - 2i - 2k - m - 1)! 



E 



Epil jym—l 



Thus, 



D(D + E) 



N _i_ h _ N ^ k m(2N -2i-2k-m-l)\ 



E 



En pZ jym—l 

DE l = D + Y} z= iE z 

DE l D rn-l = D -m-l+l + E z D^ 1 



N-i-k 

E 



m(2iV -2i-2k - m- l)\ 



- (JV_i-fc)!(JV-»-Jb-m)!^ 
It can be verified that 



m—l 



2 = 1 



N-i-k 



E k ~ 1 D (D + E) 
l 



N-i-k 



\ - m (2jV - 2t - 2A: - m - 1)! ^ / fc-1 m _ m x - 
^ (JV-i-fc)!(JV-i-fc-m)!fM ^ 

Tt=l v ' v ' 1=0 \ 2=1 



-i-k—ljjm—l 



Furthermore, 



D E k ~ l D (D + E) 



N-i-k 



DE k - 1 =D + J2 k w - = \E W 
D E z+k - 1 = D + £*t fc f 1 E v 

N-i-k 



E 



m(2N -2i-2k-m-l)\ 
\ (N — i — k)\(N — i — k — m)\ 

w rMn—l \ 



D m-l+2 + E w D m ~ l+1 + I D m ~ l+1 + I EWD 
1=0 V w=l w=l 

Another use of the expression (||) gives 

(D + Ef- 1 = Y - P(2 ^~ P ~ 3)! N , T EiDV-o 
fe(i-l)Ki-P-l)!^ 



Let now 



7 = 7(m, fc, i, N) = 
5 = 5(i,p) 



m(2N -2i-2k-m-\)\ 
(N — i — k)\(N — i — k — m)\ 

p(2i — p — 3)! 



(i - l)\{i-p- 1)! 
Applying (|2l|) , ([22]) , you can transform the equations (|l9|) and (^0|) to 



(19) 



(20) 

(21) 
(22) 
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Hence, 



DE k - 1 D(D + E) N ~ i ~ k = 

N—i—k m I k-1 z+k—1 ^ 

12 T ( £> m ~ Z+2 + ^ E w D m - l+1 + Z L> m -' +1 + / ^ E w D m ~ l 

m=l 1=0 \ w=l w=l / 

i-l p 

(D + Ef- 1 = Y12 6 EqDP ~ q 

p=l q=Q 



i— 1 p N—i—k m 

{D + E) l - 1 DE k - l D{D + E) N - i ~ k = Y^Y. 12 12 6 ^ 

p=X g =o m=l 1=0 

( E q D p ~ q+m ~ l+2 + E q D^ p ~ q ~Ep D m ~ l+1 + 

V u;=i 

z+fc-1 , \ 

Z E q D p ~ q+m - l+l + I J2 Eq D p ~ q E w D m - 1 (23) 

UJ=1 / 



In the expression fl23| ) we can change the products indicated by the horizontal braces. Using 
(|l~8|) one can obtain 

nP-iv w ( P-q~a + w-l \ r a,x? ( P-q~b + w-l\ b 

(24) 

Note that ^3) is valid only for p > q. The remaining arrangement is the substitution of 



the (24) into the equation (23). With the help of the pre-f actors 

Ax-- 



p — q — a + w — 1 
p-q-1 



a ( p — q — b + w — 1 
\ it? — 1 



the expression (|23|) turns into the form 



i— 1 p— 1 N-i—k in 

(D + E) i - 1 DE k - 1 D(D + E) N - i - k = 1212 12 12 6 ^ 

p=l q=0 m=l 1=0 

(k—1 w 
j^q £)P-q+m-l+2 _j_ j jjp-q+m-l+l _|_ _|_ 
w=l a=l 

k—lp—q z+k—l w z+k—lp—q \ 

iu=l 6=1 mj=1 a=l u;=l fe=l / 

i-l N-i-k m ( k-1 z+k+1 

+ 12 12 12 5 ^[ E p D m - l+2 + Y E p+w D m ~ l+1 + I E p D m - l+1 + I Y E p+w D r - 

p=l m=l 1=0 \ w=l w=l 



-\iri—l 



which together with the rules (|5|),(|6|) yields 



, i-lp-lN-i-k m ( , i y q , 1 y p-q+m-l+2 / i y q / i \ p~q+m-l+l 

iu=la=l ' ui=l 6=1 ' 

z+fc-1 w , -, y g+a / -i \ m— Z z+fc-lp-<2 / 1 \ o / i \ 6+m-O 

u;=l a=l x ^ ui=l 6=1 \ / \;-/ j 

+ 5 s e e/t - y + e - y + 

p=l m=l 1=0 \ w=l ' 

/ 1 \ V / 1 \ m— 2+1 2+fc+l / i \ p+mj / i \ m— n 

G) +' (i) }<*> 
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